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This paper studies a class ofZ) = n + 2(>6) dimensional solutions to Lovelock gravity that 
is described by the warped product of a two-dimensional Lorentzian metric and an n-dimensional 
Finstein space. Assuming that the angular part of the stress-energy tensor is proportional to the 
Finsteln metric, it turns out that the Weyl curvature of an Finstein space must obey two kinds of 
algebraic conditions. We present some exact solutions satisfying these conditions. We further define 
the quasilocal mass corresponding to the Misner-Sharp mass in general relativity. It is found that 
the quasilocal mass is constructed out of the Kodama flux and satisfies the unified first law and the 
monotonicity property under the dominant energy condition. Making use of the quasilocal mass, 
we show Birkhoff’s theorem and address various aspects of dynamical black holes characterized by 
trapping horizons. 


I. INTRODUCTION 

In general relativity, or more general gravitational theories admitting diffeomorphism invariance, the role played by 
the mass is quite different from that in other branches of physics. The weak equivalence principle makes it impossible 
to construct a well-defined “local” gravitational mass since it is always possible to set the local gravitational energy 
to vanish by working in a local Lorentz frame. This conceptual obstacle forced us to focus upon globally defined 
conserved quantities, e.g. the Arnowitt-Deser-Misner (ADM) mass in an isolated system [1]. However, one can 
circumvent this difficulty when the spacetime admits a spherical symmetry. In this case, the gravitational degrees of 
freedom is localized because of the absence of gravitational waves and it turns out to be extremely useful to define the 
“quasilocal” mass referring to the compact and orientable surfaces, as first demonstrated by Misner and Sharp [2]. 
A number of local geometric properties of spacetime are encoded in the Misner-Sharp mass [3], including the causal 
property of central singularities, trapped surfaces and asymptotic charges. It follows that the Misner-Sharp mass 
provides a fruitful venue for dexterously capturing the dynamical aspects of gravitational collapse. 

Inspired by recent advances of string theory, many people have tried to extend Einstein’s gravity. A covariant grav¬ 
itational theory constructed by Lovelock [4] is a natural extension of general relativity into D{> 5) dimensions. The 
most appealing and characteristic feature of Lovelock gravity inherited from general relativity is that the field equa¬ 
tions continue to remain second order, irrespective of the fact that they are accompanied by higher-order polynomials 
of curvature tensors. This traces back to the topological interpretation of each Lovelock term as the dimensionally 
continued Euler densities, allowing the understanding of Lovelock terms in the context of the BRST cohomology [5]. 
Therefore, there appear no ghost degrees of freedom at the linearized level [6, 7], and the Lovelock gravity is a 
classically well-posed gravitational theory (see [8, 9] for a recent discussion opposing this belief). Apart from this 
theoretical aesthetic beauty, the quadratic Lovelock term dubbed the “Gauss-Bonnet” term arises as a low-energy 
effective action in heterotic string theory [10-12]. This motivates us to explore quantum aspects of higher-curvature 
terms as well in light of AdS/CFT correspondence [13]. 

Since higher-curvature terms come into play where the gravitational force becomes very strong, black holes are the 
best test beds in which deviations from general relativity are capitally encoded. The complexity of Lovelock field 
equations has restricted the analysis especially to the spacetimes with a high degree of symmetry. Among other things, 
many works have focused upon the spacetime which is the warped product of two-dimensional Lorentzian spacetime 
and an n-dimensional maximally symmetric space. The simplest solution is the spherically symmetric black hole 
found by Whitt [14], as a generalization of the Schwarzschild solution in general relativity. Thermodynamics [15, 16] 
and gravitational instabilities [17-22] of this type of black hole have been intensively studied. This class of metrics 
also contains the Tolman-Bondi inhomogeneous dust spacetime [23-25] and the Vaidya-type radiating solution [26- 
29], both of which describe the gravitational collapse. The examinations of gravitational collapse have revealed that 
the global structure turns out to be quite different from that encountered in general relativity, and a peculiar type 
of massive singularity emerges in every odd dimension. In these analyses, the generalized Misner-Sharp quasilocal 
mass [30, 31] plays an essential role as in general relativity. In this sense, the Misner-Sharp quasilocal mass is more 
advantageous than the Brown-York quasilocal mass [32] constructed based upon the Hamiltonian formalism. 

In general relativity, the n-dimensional maximally symmetric space can be replaced by arbitrary Einstein spaces, 
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since their Weyl curvature fails to contribute to Einstein’s equations.^ In Lovelock gravity, on the other hand, the 
Weyl tensor appears explicitly in field equations and, therefore, the generic Einstein manifold fails to satisfy the 
vacuum field equations of the warped metric. The condition that this type of metric admit vacuum solutions imposes 
two conditions upon the Weyl curvature of the Einstein space. When one takes the Einstein space satisfying these 
conditions, the causal structures of the black hole considerably differ from those with maximally symmetric horizons, 
as argued in the Gauss-Bonnet gravity [34-36] and in the third-order Lovelock gravity [37]. This motivates our present 
attempt to explore the conditions for an Einstein horizon in general Lovelock gravity, by extending the analysis in 
[34-37] 2. 

In this paper, we generalize the previous studies [34-37] into Lovelock gravity where the spacetime consists of the 
warped product of the two-dimensional Lorentzian metric and the n-dimensional Einstein space. We find that the 
Weyl tensor of n-dimensional Einstein space must obey two conditions. We find that a dozen new Einstein spaces 
turn out to satisfy the Lovelock field equations. We extend the definition of the Misner-Sharp type quasilocal mass 
adapted to the present context. It turns out that the quasilocal mass displays some desirable physical properties 
under suitable energy conditions, provided that some conditions on the Weyl curvature and the coupling coefficients 
of Lovelock action are satisfied. Using the quasilocal mass, we further explore the properties of trapping horizons and 
their thermodynamics. 

The present paper proceeds as follows. In the next section, we give a brief review of Lovelock gravity and derive 
field equations under the setup described above. In Sec. Ill, we define a quasilocal mass as a generalization of the 
Misner-Sharp mass. We explore a number of properties of dynamical black holes defined by the trapping horizons 
using the quasilocal mass in Sec. IV. Einal remarks are described in Sec. V. In the Appendix, we give a variety of 
exact solutions for vacuum and electrovacuum cases. We follow the conventions of Wald’s textbook [42] for curvature 
tensors. 


II. SETUP 


The action of Lovelock gravity in D dimensions is [4] 


S = - d^xv^ 
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( 2 . 1 ) 


where am are real constants and we set ai = 1 and 8ttG =1. fc is given by fc = [{D — 1)/2J, where the symbol [ccj 
denotes the integer part of x.^ Sma,t is the action for the matter field, A is a cosmological constant and S denotes the 
totally antisymmetric product of Kronecker delta normalized by 
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The gravitational field equations derived from the action (2.1) reads 

Qtiv = , (2.3) 

where = —25Smat/5g^'' describes the stress tensor of the matter fields. The Lovelock tensor is given by 
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(2.4) 


obeying the Bianchi identity = 0. A notable feature of Lovelock gravity is that the equations of motion involve 

no more than the third derivative of the metric. 

In this paper, we consider the D = n + 2-dimensional spacetime for which the metric takes the cross 

product of the two-dimensional orbit spacetime {M'^,gab) and the n-dimensional Einstein space Namely, 

the local metric reads 


ds^ = gabiy)dy‘'dy^ -f P{y)j^j{x'")dx"dx^ , 


(2.5) 


^ The replacement to the Einstein space has a significant impact upon the linear instability of black holes [33]. 

^ The solutions to special cases of Lovelock gravity with a more general base space were also studied in [38—41]. 

^ Note that some literatures have employed the convention k = [L^/2j different from ours. Since the Dl2t]i term in even dimensions 
amounts to the topological invariant, it fails to contribute to the field equation. Hence, both of these conventions do not make any 
physical difference. 
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where r is the scalar on corresponding to the warp factor. Indices a, b,... run over 0,1 and i,j,... correspond to 
those of the Einstein space. The Ricci tensor of Einstein space reads RijYl] = (n — l)K'yij, where k is the 

constant normalized by k = ±1,0. We assume that (]C^,^ij) is compact with the area and that {IVP.gab) is a 
time-orientable Lorentzian manifold. The Riemann tensor of (2.5) decomposes into 

Rated ^ ^ Rated 7 Raitj 7 Rijkl ^ [b] ‘ 2 (]Dr'j 7 (^■^) 

where Da is a covariant derivative with respect to gat and (Dr)'^ = g°‘^{Dar){Dtr). The suffix “2” is attached with 
quantities of which are distinguished by those of represented by [ 7 ]. One can express the Weyl tensor 

of the Einstein space as 

^./'[7] = a/'[7]±<'- (2.7) 

When is maximally symmetric, the present setup collapses to the case analyzed in [31]. Note that the 

Einstein space is necessarily maximally symmetric in n = 3, the nontriviality arises in > 6 dimensions. 


A. Field equations 

In the following calculation, we frequently use the quantities lE(s)*j and W{s), which are defined as 
W(sY =l^i S — 0 


3Jlj2---32s-lj2s 
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W{s)ij are symmetric tensors on (/C", 7 ij). We also define 
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By a straightforward computation, the Lovelock tensor decomposes into 
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For a generic Einstein space (/C", 7 ij), W{m) are functions dependent on the coordinates and the (trace-free part 
of) symmetric tensor W{m)ij is nontrivial. In that case, the Lovelock tensor involves the convoluted coordinate 
dependence on x®, as well as the dependence on y°‘. In order to avoid these technical difficulties and make the 
discussion focused, we impose in this paper the following two conditions on (/C", 7 ij): 


77 _ /TT) 

W{mYj = —-— 6^jW{m) , 

W(m) = const. 


(2.13a) 

(2.13b) 


With these conditions, the x® dependence of drops out except for the contribution stemming from the metric 
"fij. In [34] a similar condition was imposed in the m = 2 case. On account of the dimensionally dependent Lovelock 
identities [43], the constraint (2.13a) is automatically satisfied for m > [{n + 1)/2J = k. Obviously, the conditions 
(2.13) restrict the permissible horizon topologies for static black holes. Appendix A illustrates some explicit examples 
of Einstein spaces satisfying (2.13). 

The stress tensor compatible with these assumptions, therefore, reads 


Tfj.^dxf^dx" = {fab{y) - P{y)gab)dy‘^dy'’ + r'^{y)p{y)-fij{x)dx'-dx^ , 


where 


p = —r“„ 


g^^fab = 0 . 


The stress-energy tensor obeys the conservation law 
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(2.15) 

(2.16) 


Under these settings, the Lovelock field equations are given by 
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If r is not constant, the angular part of field equation (2.19) follows from (2.17), (2.18) and (2.16). 

In the above we assumed (2.13) to simplify the system. As far as the vacuum solution is concerned, the condition 
(2.13) actually follows from the consistency with the field equations, as shown in [34] for Gauss-Bonnet gravity. 

In the following sections, we shall discuss the general properties of the metric under suitable energy conditions, 
without resorting to the exact solutions. In Appendix B, we give some exact solutions of physical interest. 


III. QUASILOCAL MASS 

In general relativity, the spacetime admitting spherical symmetry allows no freedom of gravitational radiations. 
This fact enables us to localize the gravitational energy and one can define the quasilocal mass [2] that plays an 
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important role in the analysis of dynamics [3]. When (/C", 7 ij) is a maximally symmetric space, analogous quasilocal 
quantities have been generalized to Gauss-Bonnet [30] and to Lovelock gravities [31]. These definitions have been 
further extended to the case of Einstein spaces in the Gauss-Bonnet gravity [35] and in the third-order Lovelock 
gravity [37]. Here we complete the series of research by studying the nonconstant curvature case in the full Lovelock 
gravity, which encompasses all the previous studies. 

By mimicking the quantity in Maeda’s paper [35], our proposed definition of the quasilocal mass reads 


M{y) = V: 
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1 an 
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(3.1) 


Due to the assumption (2.13b), one may view the quasilocal mass as a scalar on (M'^,gab)- It is constructed out of 
the areal radius r and its first derivative, as well as the Weyl tensor of Einstein space. When the space (/C”, 7 y) is 
maximally symmetric, the above definition reduces to the one given in [31]. It also recovers the well-defined Misner- 
Sharp mass [2] in general relativity (am >2 = 0), and its generalization in Gauss-Bonnet gravity {am >3 = 0) [30, 35]. 
Using (2.11), one can easily verify that M satisfies the variation formula 

D^M = V:r^{<3abD\-g\Dar). (3.2) 

This formula takes exactly the same form as those analyzed in previous studies [30, 31, 35]. This relation is of crucial 
importance in the following discussion. 


A. Locally conserved currents 


The physical meaning of M is less clear in the geometric definition (3.1). In this section, we shall demonstrate that 
M can be rebuilt in terms of a locally conserved energy flux. 

To proceed, let us first define the Kodama vector [44] 

AT'" =, (3.3) 

where = f-ab{(ly °‘)and tab is a volume element of {M^,gab)- This current can be viewed as a vector field 
on since AT* = 0. The Kodama vector fulfills the following crucial property: 

iL^A^ = -(Vr)2. (3.4) 

This means that the Kodama vector is timelike (spacelike) in the untrapped (trapped) region and specifies the preferred 
time direction. Using the Kodama vector, one can also define a Kodama current, 

= -gf^,K'', (3.5) 

which is again a vector field on M^. Using the Lovelock field equation (2.3), one obtains Hence 

describes an energy flux. Because of the properties K°'Dar = 0 and gabD°‘K^ = 0, one sees that these vectors are 
divergence free 

V^A'^ = 0, V^,r = 0. (3.6) 

Using DaT = —EabK^, one can easily derive the following relations, 

A:“ = -r-'^+^^Dbiy/V ^), = -r-'^yDb{M/Vy , (3.7) 


where 

VK. 

n -I- 1 


(3.8) 


is a weighted volume of /C”. It follows that vector quantities and are the Hamiltonian vector fields with the 
corresponding Hamiltonian V and M. This expression makes the divergence-free property (3.6) rather manifest. It is 
also obvious to see 


u = - y A'^u^dE , 




M = - 


(3.9) 
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where E is a (Z? — l)-dimensional hypersurface without an interior boundary and is a future pointing unit normal 
to E . This accomplishes our first aim to prove that M is a quasilocal quantity associated with the locally conserved 
energy flux. The definition of the Misner-Sharp quasilocal mass based upon the energy flux illustrates the direct 
physical relevance rather than the original geometric definition (3.1). 


B. Unified first law 

The first law of thermodynamics is one of the fundamental laws in nature. Hence, the validity of the first law 
deserves a nice criterion for the well-defined mass. We can easily check from (3.2) that the following unified first law 
[45] holds 


dM = + PdV, (3.10) 

where P has been defined in (2.15), and 

, A = Vppr^. (3.11) 

A is the weighted area of Einstein space and is related to the volume (3.9) as DaV = ADaV. As it turns out from the 
analysis of the next subsection, 'ip°‘ describes a momentum flux. Therefore, (3.10) represents the physical circumstance 
that the energy balance is compensated by the work term PdV and the energy inflow provided by ip°‘. The unified 
1st law allows us to interpret M as an energy contained in the closed surface enclosed by the geometric radius r. 


C. BirkhofF’s theorem 

The theorem of Birkhoff plays a significant role when one analyzes the gravitational collapse of a spherical body 
in general relativity. One of the characteristic features of Lovelock gravity is that BirkhofF’s theorem (and modified 
versions thereof) continues to hold, as discussed in [30, 31, 46-48]. 

Let us consider the matter fields with Tab = 0. Suppose that the first line of equation (2.17) is nonvanishing.This 
implies that 

0 = DaDbV - ^D'^rgab ■ (3.12) 

Assume that DaV does not vanish. Then, DaV describes a conformal Killing field on which implies 

DaKb = ^D‘^reab. (3.13) 

It follows that is a Killing vector on M^. Thanks to the property = DaKb{dy°')fi{dy^),y, describes a 

hypersurface-orthogonal Killing vector on M. , 

X[^V.iLp]=0, V(^iL,)=0. (3.14) 

On account of the property (3.4), this means that the spacetime is static in the untrapped region. Hence, Tab = 0 
provides a sufficient condition for the validity of staticity (in the untrapped region). This also justifies the physical 
interpretation of as a flux current, since it vanishes in static spacetimes. 

In order to obtain the metric explicitly, let us concentrate on the vacuum spacetime (P = p = 0) in what follows. 

Then, the variation formula (3.2) [or the unified first law (3.10)] immediately gives M = g = const. Let us introduce 

the coordinate thy K = d/dt and denote the norm of by —/, i.e, Under the condition (Dr)'^ ^ 0 

one can use r as the coordinate on conjugate to t, and the general solution reads 

dr^ 

ds2 = -f{r)dt^ + — , (3.15) 

f[r) 


We shall not discuss in this paper the case where the first line of (2.17) identically vanishes. For these artfully chosen values of am with 
a given Einstein space, there appears a solution for which the metric g^b is undetermined. For this class of metrics, we refer the reader 
to Refs. [30, 31, 46]. 
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where /(r) satisfies 

k 


^^=v: 


E 

,m—l 
k m — 1 


1 am r 


,n—2m+l 


EE 2 

m—1 l—O 


2m+i ^ ^ _ 2m + 1 
„n—2?71+1 


W{m) — 


j.n+1 

n + 1 
21 


-A 


l — r 


i + i 


( z ^ (2m- 2) +p)j (k - f{r)y+^W{m-l - 1) 


(3.16) 


This is a simple Lovelock generalization of the Dotti-Gleiser solution [34] for Einstein-Gauss-Bonnet gravity. 

The case r = ro = const, also solves (3.12). In this case, (2.19) gives that is constant, thus {M^,gab) is a 
spacetime of constant curvature. The metric can therefore be written as 


ds^ = -(1 - Aa;^)dt^ + --^ + rljijdx'dx^ , 

1 — Xx^ 


(3.17) 


where tq and A satisfy the following relations 
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(3.17) describes a Nariai-type metric (A)dS 2 x /C". 


m-l , . / 2 i + l \ 

EG’,‘ n (n — 2m + 1 + p) j K^W{m — l — l) 

1=0 ^ ^ \p =0 / 

T ) ~ k’‘W{ m — 1). 


2 (n — 2m + 21) f m 
n{n — 2m) 
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D. Physical properties of quasilocal mass 

In order to analyze the dynamics of the warped spacetime (2.5), it is advantageous to work in the double null 
coordinates 


ds^ = —2e ^^^’'^'^dudv + r'^{u,v)'yijdx'^dx\ (3.20) 

where the orientation is fixed to be euv > 0. Then we can see that the variation formula (3.2) may be cast into 

duM = -v:efr'-+^ (T„,0_ - r„„0+) , (3.21a) 

n 

d,M = ilG=e^r"+i (T„,0+ - T,,9_) , (3.21b) 

n 

where 9± describe the expansion rate for the null directions 

(3.22) 

r r 

The variation formula (3.21) does not involve the Lovelock coupling coefficients nor the information on the Weyl tensor 
of the Einstein space explicitly. This fact is advantageous for discussing the monotonicity property of the quasilocal 
mass as described below. 

We fix the spacetime orientation by declaring that the future-directed null vector d/dv (resp. d/du) is outgoing 
(reps, ingoing). Namely, 0+ > 0 and 0_ < 0 hold on an untrapped surface. Remark that each value of 0± is not an 
invariant quantity by virtue of the remaining freedom of rescaling u —> U{u),v —> V{v). Instead, e^6+9- enjoys an 
invariant physical meaning characterizing the trapping nature. 
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In order to extract the physically reasonable results, we impose the energy conditions on the matter fields. The 
null energy condition for the matter field implies 

Tun >0, T™ > 0 , (3.23) 

while the dominant energy condition for the matter field implies 

Tun > 0, T™ >0, T,, > 0. (3.24) 

Note that only the information of radial directions is encoded on these inequalities. 

Let us now establish that our quasilocal mass exhibits a monotonicity property, which is desirable for M as a 
physically reasonable mass function. 

Proposition 1 (Monotonicity) If the dominant energy condition holds, the quasilocal mass is nondecreasing along 
outgoing null or spacelike directions on an untrapped surface. 

The proof follows immediately from (3.21). □ 

Let us next move on the the positivity claim. To this end, let us first define the regular center. The central point 
is said to be regular center if 

K - {Drf - Cr^ (3.25) 


holds where C is a nonvanishing constant. 


Proposition 2 (Positivity-I) If the dominant energy condition holds and the spacetime has a regular center which 
is surrounded by untrapped surfaces, quasilocal mass is non-negative. 


Suppose that W(k) is nonvanishing. Around the regular center, quasilocal mass behaves as 


M ~ v: 


1 flfc 


„n—2fc+l 


2 ^=+! k (n - 2fc + 1) 
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1 Ofc r" 
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Combining with the monotonicity property and (3.27), the dominant energy condition requires 

akW{k) > 0 . 


(3.26) 

(3.27) 

(3.28) 

(3.29) 


This proves the positivity of quasilocal mass around the regular center. The monotonicity property establishes the 
claim, as we desired. □ 

It is worth noting that for the k = 1 case, the regular center is always surrounded by untrapped surfaces, while this 
is not the case for k = —1. Remark also that the Misner-Sharp mass behaves M oc around the regular center 
for the case with /C" being the maximally symmetric space, whereas M oc 7 -"’- 2 fe+i for the present case. 

Next, let us consider the case in which the spatial hypersurface admits a marginal surface as its inner boundary. 
On the marginal surface we have (Dr)^ = 0; hence, the following version of the positivity holds. 

Proposition 3 (Positivity-II) Suppose the dominant energy condition and A < 0. If the spacelike hypersurface 
admits a marginal surface as its inner boundary, then the quasilocal mass admits a positive lower bound, provided that 
the Lovelock coefficients and Weyl tensor satisfy the following conditions for all m, 



]^(n - {2m - 2) +p) 
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> 0 . 


(3.30) 
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This is clear from the monotonicity and the definition of quasilocal mass. □ 

We obtained a condition (3.30) under which the positivity of the mass holds. Inspired by string theory, we may 
physically fix some of the Lovelock coefficients. However, it appears that the sign of Weyl tensor of the Einstein space 
is not controllable. It would be better if we have a clearer physical and mathematical meaning of (3.30). We leave 
this for future investigations. 

To conclude this section, let us make a brief comment on the asymptotic behavior of the quasilocal mass. If 
the Einstein space is a round sphere, the metric falls into the standard definition of asymptotic flatness, and it 
would be meaningful if the asymptotic value of the quasilocal mass converges to the ADM mass, as argued in 
[30, 31]. If the Einstein space is not the maximally symmetric space, the metric exhibits a slow falloff and it does not 
allow asymptotically flat/AdS solutions in the standard sense. For this reason, we shall not attempt to discuss the 
asymptotics for the quasilocal mass. 


IV. TRAPPING HORIZONS 


The concept of event horizon is not of practical use because the identification of its locus requires the knowledge 
of the evolution of Einstein’s equations into the entire future. A more convenient manner to characterize locally the 
strong gravity is the trapping horizon, which was originally proposed by Hayward [49]. In this section, we address 
some properties of trapping horizons in the present settings. 

The trapping horizons are the n + 1-dimensional hypersurface foliated by n-dimensional marginal surfaces on which 
0+9- = 0 is satisfied. Set 0+ = 0 on the marginal surface in what follows. Then the marginal surface is said to be 
future for 0- < 0, past for 9- > 0, outer for du9+ < 0 and inner for du9+ > 0. By definition, the notion of trapping 
horizons is quasilocal and does not make any references to the asymptotic structure. One may deduce intuitively that 
the future-outer trapping horizons are of the most relevance for a local description of dynamical black holes, since 
inside the trapping horizon both of the outgoing and ingoing rays are converging. In the following discussion, we shall 
be mainly interested in (future-)outer trapping horizons. According to the proposition 12.2.4 of [42], trapped regions 
cannot be causally connected to null infinity, provided the null convergence condition and the cosmic censorship are 
valid. Therefore, the existence of the trapped regions implies the event horizon in a physically reasonable condition. 

The properties of trapping horizons have been analyzed in detail for the Gauss-Bonnet gravity [50] and for the 
Lovelock gravity [31] with the maximally symmetric horizons. In order for the trapping horizons to inherit properties 
in general relativity, we have to assume a certain inequality involving the Lovelock coefficients and the Weyl tensor 
for the Einstein space. 

The next proposition specifies the causal character of the trapping horizon. The proof is the same as in Ref. [50]. 


Proposition 4 (Signature law) Under the null energy condition the outer trapping horizon is nontimelike, provided 
that the following condition on Weyl tensor and the Lovelock coefficients holds for all m, 


2 i+i - 2) + p) j k'W( m - 1 - 1) 

.1=0 ^ ^ \p=0 / 


> 0 . 


(4.1) 


Let f be a generator of the outer trapping horizon at which 9+ = 0 and i9„0+ < 0. Since the trapping 

horizon is foliated by marginal surfaces, we have 


£^9+ = Cd^9+ + Cdu9+ = 0. (4.2) 

Evaluating the (v,v) component of (2.17) at the trapping horizon we get 

^ n ol — m + l / ™ _ -1 \ / \ 

Tvv = -J2Y1 ^ 2^-2 ( I ) [Y[in-2m + 2+p)\ K^Wim-l-l)d,0+. (4.3) 

m=l 1=0 ^ ' \p=0 } 

The null energy condition and the inequality (4.1) thus assures d„9+ < 0. Hence (4.2) implies that < 0 is 
satisfied. If the trapping horizon is timelike, this inequality is not satisfied. We therefore arrive at the claim. □ 
The most interesting property of the event horizon of a black hole is the area increasing theorem (Proposition of 
12.2.6 of [42]). It turns out that a similar property holds for the trapping horizon. 


Proposition 5 (Area law) Under the null energy condition and the conditions in Proposition 4, the area of outer 
trapping horizon, A{r) = Vffr'^, increases along its generator. 
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The proof directly follows from 

C^A = nr^-^V: {^duT + Cd.r) 

= r'^V^O.C > 0 , (4.4) 

where we have used > 0 for the nonspacelike (spacelike) trapping horizon to be future-pointing (outgoing), hence 
< 0 from the signature law. This completes the proof. □ 


A. Dynamics of trapping horizon 


In general relativity, the trapping horizons display laws analogous to ordinary black-hole thermodynamics even in a 
dynamical circumstance [49]. Since the unified first law (3.10) represents the energy balance, it can be used to deduce 
the thermodynamic first law for a trapping horizon. One can recast (3.10) into 


Aljja =DaM 


k m — 1 




m—1 l—O 


{l + l)r 


X f (^ ~ (2m — 2) -I- p) j (k — {Dr)^yW{m — 1 — 1) 

\p=0 J 




m—1 


2m ^ j,2m " “ 


(4.5) 


Let be a generator of the trapping horizon. Hence, along the trapping horizon (Dr)^ = 0, we get 


where we have defined 


k m — 1 


_m—l l—O 


„n —2m+2 


n — 2m + 2 


I n ~ ~ 2) -I- p)K'‘W{m — 1 — 1) 


kth = -^D r 


(4.6) 


(4.7) 


One can interpret (4.7) as a surface gravity of a trapping horizon, since it fulfills [45] 


= KTHi^h , (4.8) 

where the equality is evaluated on the trapping horizon. Note that this equation resembles the equation defining the 
surface gravity of a Killing horizon [42]. It deserves to emphasize that the surface gravity is not constant over the 
trapping horizon, as can be inferred from the Vaidya-type radiating solution (see Appendix B.). 

The unified first law reads 5^M = Ai^tj} + PS^V, hence Aid'll) term should be identified as TS^S term. Assuming 
that the temperature is related to kth by T = KTH/(27r), we can identify the entropy of a trapping horizon as 


S = 2nVJ^ 


k m — 1 

EE 

m—1 l—O 


n-2m-\-2 


CLm^ 


l — m-\-l 


n — 2m + 2 


^ ~ ~ + 7 *)^ (to — 1 — Z) 


(4.9) 


Eq. (4.6) also justifies the physical interpretation of as flux current, since the change of trapping horizon entropy is 
responsible for the flux through the horizon. In the general relativistic case, the entropy is proportional to the area of 
the trapping horizon. The Lovelock black holes therefore admit a correction arising from higher-curvature terms [15]. 
The highest term Z = m — 1 is also present for the maximally symmetric horizons, while the other terms represent the 
contribution coming from the Einstein horizon. 

Now the entropy of a trapping horizon is obtained, we move to prove the entropy increasing law. This corresponds 
to the second law of black hole dynamics. 


Proposition 6 (Entropy law) Under the null energy condition and the conditions in Proposition 4, the entropy of 
outer trapping horizon increases along its generator. 
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The variation of the entropy along the generator gives 

k m — 1 


C^S = 2itV: 


EE 2 /-m+l^n- 2 m +2 


f n“ 2) + p)'j k'-W{ m - I - 1) 


_m—l l—O ^ ' \p—0 

The proof follows immediately from the same argument as the area theorem. 


9_C 


(4.10) 


□ 


B. Wald’s entropy 


In the previous subsection, we derived the entropy of a trapping horizon by requiring the first law of thermodynamics 
for the trapping horizon. Here we reproduce it by Wald’s prescription for the Killing horizons [51, 52]. 

Suppose that the metric admits a nondegenerate, bifurcate Killing horizon r = rt with a bifurcation surface B. 

Sw = -27r J e^.epAdC , (4.11) 

where is the binormal to B given by (3.3), 


Sw = -27r / ^ 


m—1 

k 


= -27r 


= 27rKf 


1 n D 1^12^2 Z? 2^2m-l2^2iz 

-L / ^^I^l2^2 -^'2m —l2^2m ' ^^P2m — lP2 

2"* m 
1 


dR 


pLVpX 


E 2 "’ {^^bdvli 


M4---M2m-lM2m p ^'3^'4 p ^2Tn- 

2/4...I/2m-l2^2m '^M3AA4 ' * ’ '^Ai2m-1 M2m 


m=l 
fc m— 1 




^p\dV^ 




.m=l 


n—2m+2 


n — 2 to + 2 


f™ n ~ ~ 2 ) + p)k^IK( m — 1 — £) 

d p=0 


(4.12) 


Therefore, we can see that the entropy we defined from the quasilocal mass coincides with Wald’s entropy. 

We have derived the expressions of entropy for the Killing horizons and see that it coincides with the stationary 
limit of trapping horizons. One can alternatively utilize the Kodama vector instead of the generator of the Killing 
horizon to directly derive for the trapping horizon as demonstrated in [53]. 


V. FINAL REMARKS 

In this paper we explored various properties of the spacetimes which are the warped product of a two-dimensional 
Lorentzian spacetime and an n-dimensional Einstein space. Assuming the form of the stress-energy tensor to be (2.14), 
we revealed that the Weyl curvature of the Einstein space must obey certain conditions (2.13). This assumption comes 
not only from the simplification, but also from the requirement that the metric (2.5) admits a vacuum solution. Some 
nontrivial examples are given in Appendix A. We found that all the isotropy irreducible spaces fulfill this property. 
Our study enlarges considerably the solution space of Lovelock gravity. 

One immediate conclusion for replacing the n-dimensional maximally symmetric subspace by the Einstein space is 
that the metric shows the fall-off behaviors different from the standard one. This means that the M — Vh diagram for 
the static black hole is much more complicated than [14]. A possible future work in this direction is to examine the 
P — V criticality of a black hole with Einstein horizons and to expose the thermodynamic phase structure. 

We then proceeded to define a quasilocal mass and explored its physical properties, by extending the previous 
works [30, 31]. The rederivation of the quasi-local mass in terms of the Kodama flux is desirable for the physical 
interpretation of the quasilocal mass. Up to the certain conditions among the Lovelock coefficients and the Weyl 
curvature of the Einstein space, it turns out that the quasilocal mass shares the same behavior as that in general 
relativity [3]. This implies that the Misner-Sharp-type quasilocal mass continues to be useful also in Lovelock gravity 
and can be utilized to obtain a coherent picture of spacetime dynamics as exemplified by gravitational collapse. We 
hope to come back to the point for the deeper mathematical and physical understanding of the conditions (3.30), 

(4.1). 

Our formulation of Lovelock solutions with the warped n-dimensional Einstein space is very robust and has plenty of 
potential applications. We expect that the geometrodynamics approach to Hamiltonian formulation of Lovelock black 
holes [54] can be extended to the case with Einstein horizons. It is also interesting to consider the effects of nonlinear 
Maxwell field [55] and higher-rank p-form fields [56], which would display an intriguing thermodynamic phase by the 
interplay with Weyl curvatures of Einstein spaces. One may also explore the generalization of C functions [57] and 
the maximal entropy principle [58] into the present context. 
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Appendix A: Examples of Einstein spaces 

We give some examples of Einstein spaces satisfying (2.13). The condition (2.13b) is trivially met for the homoge¬ 
neous spaces G/H since they admit a frame in which the curvature tensors take constant values. Remark that not 
all the homogeneous spaces fulfill (2.13a). For instance, we find that the homogeneous Sasaki-Einstein space 
described by the coset SU(2) x SU(2)/U(1) [59] fails to satisfy the m = 2 condition of (2.13a). It has been also shown 
in [34] that the family of homogeneous Bohm metrics does not satisfy (2.13b). The following examples are either the 
products of maximally symmetric spaces or the symmetric spaces. 


1. Products of maximally symmetric spaces 

a. X 

Let us consider the space consisting of the two products of maximally symmetric spaces x K‘^ (p + Q = n), 
where denotes the p-dimensional maximally symmetric space with a sectional curvature k = 0, ±1. Then the 
n-dimensional Einstein metric reads 




P- 1 
P + 9-1 


dfi2(A:p) 


9-1 
P + 9-1 




(Al) 


Decomposing indices into for ffP and {/, J,...} for the nonvanishing components of the Weyl 

tensor are 


ri CD 
CAB = 


Kq 

P- 1 


^CD 
°AB • 


Cai^^ = , 




Kp 

7^ 


:CD 
°AB ■ 


(A2) 


In order to see whether the constraints (2.13a) are satisfied, it is enough to check that IF(m)® is proportional to J® 
for all m. The constraint for m = I is trivial, whereas the m = 2 imposes the condition obtained in [34] and requires 
p = q. For arbitrary order, the left-hand side of (2.13a) reads 


w(m)] 


m m — s 


y2m— 


s—0 t—0 


m \ m — s 


s-\-t 


(_l)™-(s+‘) 


p - 1 / [p - (to -I- s - f) - I] [p - (/ - s -I- f)] 




(p - (to -I- s - f) - 1 -I- a) j {p - {m - s + t) + P) \ . 


(A3) 


Q;=0 


3=0 


This proves that the constraints (2.13a) are satisfied at arbitrary order of to. 


b. X ■■■ X K” 

The gth products of maximally symmetric spaces also satisfy the condition (2.13). The metric in this space is 
given by 


CT=1 / 


(A4) 


The Weyl tensor is computed to give 


a 


l(a-) j{cr) 


fcMjM _ ,^P(9- l)^feM/M 




p-1 


/(CT)j(cr) , 






for <J ^ p, 


(A5) 
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where • ■ • } are indices for After lengthy calculations, the left-hand side of (2.13a) is 


Wirnf. =k^51 


pq 


pq — 2m 


E 


! / 2 p(g- 


ni<C<J7<g “(?)(')) ■ \ P ^ 




p (a(7)(7) + S/3=l ®(/3)(7)) 1 

1 


X n 

Q!=l,a 7^7 


P {o^{oL){a.) + S/3^1'^(a)(/9)) 


n 

X~I =0 
“(c)(c)-|-Ep = l “(c)0) 

n 

Xa=0 


p ^a(7)(7) + a(/3)(7) J 1 + X 7 

P— ( 0(a)(a) + E/®(“)(/5) j 


(A 6 ) 


where a(cr)(p) (1 < ct, p < g) are numbers satisfying J 2 i<( 7 <p<q^i<y){p) ~ easily see that the right-hand 

side of (A 6 ) is independent of any choice of 7 . This proves that the constraints (2.13a) are satisfied at arbitrary order 
of m. 


2. Isotropy irreducible spaces 

We discuss here a class of homogeneous spaces satisfying (2.13). Let us start by introducing some mathematical 
definitions (see [60] for details). Let G be the semisimple Lie group. A manifold is said to be homogeneous when 
G acts transitively. Letting H denote the isotropy subgroup of G, the homogeneous space is described by a coset 
space G/H. The homogeneous space is said to be symmetric if there exists an involutive automorphism a such that 
a{h) = h, (j{k) = —k, where h is the Lie algebra of H and k is its complement in the Lie algebra g of G; i.e. we have a 
Cartan decomposition, g = h(B k. The homogeneous space is isotropy irreducible if the linear isotropy representation 
of H is irreducible and the symmetric isotropy irreducible space is simply called an irreducible symmetric space. 

As we commented earlier, the first condition (2.13a) is always satished for any homogeneous spaces G/H, whereas 
not all the homogeneous spaces meet the second condition (2.13b). Let us here remind the theorem that any symmetric 
G-invariant rank-two tensors are unique up to the multiplication constant in the isotropy irreducible spaces [61]. It 
follows that the symmetric tensor as well as the Ricci tensor, constructed out of the G-invariant metric must 

be proportional to the G-invariant metric in the isotropy irreducible spaces. Therefore, we have proven that all the 
isotropy irreducible G-homogeneous spaces are Einstein and satisfy the required conditions for the Weyl tensor (2.13). 
In what follows, we shall provide some examples by focusing on the Riemannian irreducible symmetric spaces® for 
which the explicit metrics are known. Of course, all the G-homogeneous spaces listed in Tables 7.102-7.107 of Besse’s 
textbook [60] provide desired examples. 


a. Complex projective space CP^ 


The complex projective space is a Kahler-Einstein symmetric space 
coordinates, the Kahler potential is given by 


N 


K = i^\nil + J2 




= SU(iV -k 1)/SU(77). In the complex 

(A7) 


7=1 


The standard Fubini-Study metric ds^ = 277 jdz^d 2 :'^ with jjj = d^K/dz^dz'^ reads 


ds^ = 


2^ 


1 + kl 


dz^dz^ - 




(AS) 


1 + I^P 

where I = 1,N = n/2. The curvature tensors in the complex basis are given by 

R'Lb] =-2r^6\jS^L), R^jb] = {N + l)e-^S^j, Rb]=2N{N+l)rb (A9) 


^ As far as the noncompact case is concerned, this exhausts all the possibilities (Proposition 7.46 of [60]). 
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Hence, = (n + 2)/[2(n — 1)] gives the unit normalization. In the real basis, we have 


a/'[ 7 ] 





-\-Ji-jJ + 


, /.[fc 7-^] 

Ul U j 


(AlO) 


where J is the complex structure satisfying Ji^Jk^ = —5i^ and Hence, it is almost obvious to see that the 

symmetric tensor W(jn)ij = W{'m)ji is proportional to ^ij for each to, although we omitted the detailed computations. 

One can also consider the Bergmann space, which is the negative curvature version of CP^; namely, it is described 
by a coset SU(7V — 1,1)/S(U(A^ — 1) x U(l)). This metric arises as nonlinear sigma models of vector and hyper 
multiplets in the framework of A/” = 2 supergravity [62]. The Kahler potential is given by 


N 


K = 




1=1 


(All) 


Since the curvature tensors of the Bergmann space are given by the Wick rotation ^ ^ ii oi (A9), it is obvious that 
the Bergmann space also satisfies (2.13), as the complex projective space. 


b. Quadratic surface 0^"^(C) = SO(A + 2)/SO(A) x SO(2) 


The quadratic surface is the compact space of positive curvature in n = 2N dimensions. The noncompact version 
S0(1V, 2)/SO(7V) X SO(2) appears as the nonlinear sigma model for A/" = 2 supergravity [63]. The Kahler potential is 
given by 


N N 

a: = £2 log I 1 + 2 ^ z^z^ + ^(^-^ 2^)2 

I I,J 


The curvature tensors in the complex basis are 

R^j^l = -r\5^j5^L + S^lS^j - , 

giving = n/[4(n — 1)]. It is then straightforward to show explicitly that the condition (2.13) holds. 


(A12) 


(A13) 


C. G2(2)/S0(4) 


This is a negative curvature, noncompact space and describes a nonlinear sigma model arising from the U(l)^ 
reduction of five-dimensional minimal supergravity. Setting ^ > In / and / —>■ —/ [this sign flip ensures H = SO(4) 
instead of 0(2,2)] in Ref. [64], the metric is given by ds^ = where 


el = 


d/ 

V2f 


e^ = 


V2f 


(dx + uidui-I-n2dM2), e'i = W-d(/, 


= \ l -e'^dK, 


e® =—=e , e® = —^=e^‘^^^[dui + K^d'Ci — \/3K{du2 — Kd'(; 2 )], 

v 2 / v 2 / 


(A14) 


e^ = 


1 


V 27 


e + VSKdvi), 


1 


e = 


vv 


e^^^{du2 — 2Kdv2 — v^K^dui). 


Here eight variables {f,x,4>,K,vi,ui,V2,U2) span the coordinate basis. The Ricci tensor is Rij[x\ = hence, 

= d/1 gives the correct normalization. From the general argument given at the beginning of this appendix, this space 
satisfies (2.13), which can also be checked explicitly (the quartic term is trivial due to the dimensionally dependent 
identities [43]). One can also confirm that the curvature tensors take constant values in the above orthogonal basis. 


d. SL(A,R)/SO(A) 


This is the negative curvature symmetric space of dimension n = -|- 2){N — 1). The coset arises as a nonlinear 

sigma model when one dimensionally reduces the H-dimensional Einstein’s gravity on an A^-dimensional torus [65]. 
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The metric reads 


f- 

7 y dxMa:'^ = — —Tr(dMdM ^), (A-15) 

where M is an x iV symmetric matrix with detM = 1. The Ricci tensor is given by = —{N/£‘^)jij, hence 

£2 = 2N/{N‘^ + N — 4:). A convenient parametrization of the matrix M is M = U]^UdUr where Ud is the diagonal 
imimodular matrix and Ur is the upper triangle matrix with the diagonal entries equal to unity. To the best of our 
knowledge, the Riemann curvature of this space does not seem to have a simple expression for general N. In spite of 
this, the group-theoretical reason ensures (2.13), as can be checked for small N explicitly. 


Appendix B: Inclusion of matter fields 

We present in this appendix some solutions with matter sources. These solutions are instructive to give some 
physical insights for readers. 


1. Maxwell field 


Let be a Faraday two-form satisfying 


di^ = 0, d*A = 0. 


with a stress-energy tensor 


T^. = 2{ - -g^^Fp^FP^ ) . 


Due to the restriction Tai = 0, the permissible form of the two-form reads 

F = ^Eeabdy^- A dy’’ + ^F,jdx^ A dx^ . 
From the field equations, we get E = E{y), Fij = Fij{x^) with 

Q 


E=—, dF = 0 , 


d F = 0 , 


(Bl) 


(B2) 


(B3) 


(B4) 


where Qe is a constant. Hence, F is a harmonic two-form on /C”. The nontrivial F exists provided the second Betti 
number of /C" is nonvanishing. Moreover, the stress-tensor of the electromagnetic field must fall into the form (2.14), 
which yields the additional restrictions 

FikFj^ oc , Fijpd = const. (B5) 


The existence of nontrivial F satisfying these conditions also puts a restriction on the Einstein manifold. For instance, 
the condition (B5) is satisfied if /C” is even dimensional and admits an almost complex structure J with Fij oc Jij. 
Note that in odd dimensions we have detFy = 0; hence, the first condition of (B5) implies Fij = 0 provided /C" is 
Euclidean (this can be shown by taking the determinant of 'y'^^FikFji = n~^')ijFkiF^^. See e.g. [66] for a similar 
analysis). 

Denoting FikFj^ = Q^'^ij where Qm is a constant, we have 


Tab 


Ql 


2r4 




T — 


Ql (4 - n)Ql 


+ 


2r2 


7ij 


(B6) 


Since Tab = 0 is fulfilled, it follows that Birkhoff’s theorem holds for the present field configurations. The solution of 
(Dr)^ > 0 is given by (3.15) with the left-hand side of (3.16) replaced by 


Remark that we have Qm 
at infinity. 


y ^ fi + V^ 


Ql 


nQl 


(1 — n)'i 


.n—1 


2{n — 3)r 


0 for odd n. One sees that the nonvanishing Q 


3w) ■ (B7) 

m gives rise to the slow falloff of the metric 
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2. Vaidya solution 

Another interesting solution of physical interest is the Vaidya-type radiating solution [26-29]. Let us consider the 
null dust fluid as the stress-energy tensor 


, (B8) 

where is a null vector and p is the energy density of the null dust. Here, let us introduce the outgoing null 
coordinates: 

dsf = —f(v, r)(iv^ + 2dndr , . (B9) 


It follows that P = 0, Tab = P^avVbV] hence, the unified first law implies AM = pAAv. This means that M is a 
function of v only and the metric function f(v,r) satisfies 


M{v) =v: 


E 

^m—l 
k m — 1 


1 


Ciri 


„n—2m+l 


+ E 

m—1 l—O 


2m+l _ 

„n—2m+l 




j-n+l 

n -I- 1 
21 


-A 


a„ 


l+l 


- 2) + p) j {k - /)'+^IT(m -1-1) 


p{v, r) is given by 


1 dM{v) 
dv 


(BIO) 


(Bll) 


The null energy condition requires that M(v) is an increasing function of v. 
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